Headgroup dimerization in methanethiol monolayers on the Au(lll) surface: a 

density functional theory study 
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A long-standing controversy related to the dimer pattern formed by S atoms in methanethiol 
(CH3SH) on the Au(lll) surface has been resolved using density functional theory. For the first 
time, dimerization of methanethiol adsorbates on the Au(lll) surface is established by computa- 
tional modeling. For methylthiolate (CH3S), it is shown that the S atoms do not dimerize at 
high coverage but reveal a dimer pattern at intermediate coverage. Molecular dynamics simulation 
at high coverage demonstrates that the observed dialkyl disulfide species are formed during the 
desorption process, and thus are not attached to the surface. 
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PACS numbers: 61.46.-w, 36.40.Cg, 68.43.Bc 
I. INTRODUCTION 

Recently, much attention has focused on the proper- 
ties of self- assembled monolayers (SAM), particularly of 
alkanethiol on gold surfaces see Ref. [HJ for a 

review. This high level of interest may be ascribed to 
their relevance to wetting phenomena, tribology, chem- 
ical and biological sensing, optics and nanotechnology. 
Despite the apparent simplicity of these systems, their 
observation in various experiments has led to contradic- 
tory results, in particular related to the question if the 
S-H hcadgroups of alkanethiol organize in dimers on the 
Au(lll) surface or not. 

In the early phase of SAM related studies, it was as- 
sumed that alkanethiol adsorbs at threefold coordinated 



hollow sites on the Au(lll) surface [15J. This picture has 
been first challenged by grazing incidence X-ray diffrac- 
tion (GIXRD) studies [l|,l2|], where a dimer pattern was 
proposed. High-resolution electron-energy-loss spectra 
supported this dimer model 0]. Scanned-energy and 
scanned-angle photoelectron diffraction experiments on 
alkanethiolate, in contrast, did not yield any dimeriza- 
tion [H, Q . Based on temperature-programmed desorp- 
tion, Auger electron spectroscopy, and scanning tunnel- 
ing microscopy (STM) studies, methanethiol was found 
to form dimers [lfj. Recently, Torrelles et al. observed 
by use of GIXRD, STM and electrochemical techniques 
that alkanethiolate does not dimerize [HI]. 

The present contribution addresses these controver- 
sial observations by simulation from first principles. A 
consistent model is proposed that agrees with all avail- 
able experimental data related to thiol adsorption on 
the Au(lll) surface. Until recently, it was assumed 
that the S-H bond of methanethiol is cleaved once it 
is adsorbed to the Auflll) surface, such that it be- 
comes methylthiolate [J, [lfj. Thus, methanethiol and 
methylthiolate would exhibit the same adsorption pat- 
tern once deposited on Au(lll). Recent results, how- 
ever, show that methanethiol stays intact when in con- 
tact with the defect-free Au(lll) surface [ll], [T^], and 



turns into methylthiolate on Au(lll) only in the pres- 
ence of vacancies. To examine if methanethiol dimerizes 
on Au(lll), it is crucial to study the methanethiol and 
methylthiolate dimer patterns separately. If these two 
species display different adsorption structures, the long- 
standing controversy related to their dimerization behav- 
ior might be resolved. 

Guided by this motivation, we study in the 
present contribution the adsorption of methanethiol and 
methylthiolate on the Au(lll) surface by use of density 
functional theory [l6l ]. First, we will present adsorption 
energies and geometries for methanethiol on the regular 
Au(lll) at 0.5 ML and 1.0 ML, and demonstrate, for the 
first time, that methanethiol reorganization on the clean 
Au(lll) surface gives rise to the formation of a dimer 
pattern. Subsequently, we will show that methylthiolate 
on the regular Au(lll) surface exhibits a coverage de- 
pendent dimerization behavior. Further, we address the 
possible impact of surface defects and perform molecu- 
lar dynamics (MD) simulations to inspect the effect of 
temperature on the dimer formation. 



II. COMPUTATIONAL METHOD 

We employ the VASP code [lj| which involves the pro- 
jector augmented wave potential Il71 . | X 811 and the PW91 
generalized gradient approximation [19j|. The wave func- 
tions are expanded in a plane wave basis with an energy 
cutoff of 400 eV. The Brillouin zone integration is per- 
formed by use of the Monkhorst-Pack scheme [2(j. We 
utilized a 3 x 3 x 1 k point mesh for the geometry op- 
timization. A (3 x 2y3) superlattice (see Fig. [1^) is 
employed as an Au(lll) supercell corresponding to 
12 Au atoms per layer. The Au atoms in the top three 
atomic layers are allowed to relax, while those in the bot- 
tom layer are fixed to simulate bulk- like termination [2l| . 
The vacuum region comprises seven atomic layers, which 
exceeds substantially the extension of the methanethiol 
molecule. To examine the accuracy of our approach we 
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TABLE I: The geometries and adsorption energies for the 
studied methanethiol configurations on the regular Au(lll) 
surface at 0.5 ML and 1.0 ML. The entries ds-s, ds-Au (in 
A) and E a d s (in eV/per methanethiol) refer to the S-S bond 
length, the shortest Au-S bond length and the adsorption en- 
ergy, respectively. The superscripts I and E stand for the ini- 
tial and equilibrium structures, respectively. The maximum 
adsorption energy is underlined. 





FIG. 1: (Color online) (a) (3 x 2\/3) supercell. (b) top+hcp. 
(c) hcp-hcp. (d) top+hcp'. (e) top+fcc+3b. (f) hcp-hcp* 
with a vacancy, (g) top+top. (h) top+fcc+3b'. 



increased the energy cutoff to 500 eV and the number of k 
points to 8 x 8 x 1. In both cases, the difference amounted 
to less than 2%. We further calculated the gold lattice 
constant, and found it to agree with the experimental 
value [H| within 2.1%. 



III. RESULTS AND DISCUSSION 

A. Methanethiol molecules deposited on the 
regular Au(lll) surface 

In what follows we present our results, beginning with 
a discussion of both the geometries and adsorption ener- 
gies for the optimized configurations of methanethiol on 
the regular Au(lll) surface at 0.5 ML (Fig. [lg, Fig. [Tfr) 
and 1.0 ML (Fig. [Tp - Fig. [if) coverage, as displayed 
in Table HI At 1.0 ML, we arrange four methanethiol 
molecules in the (3 x 2-\/3) supercell according to the ex- 
perimentally detected structure [H, [24| . The molecules 
labeled 1(3) and 2(4) are symmetry equivalent, see Fig. 
[T^. The nomenclature top+hcp in TableUJrefers to an ini- 
tial geometry where the S atom of adsorbate 1 is placed 
on an exact top site, while that of adsorbate 3 is just in 
an fee region so that the distance between the S atoms 
of 1 and 3 is close to 2.3 A (see Fig. [Tp), and analo- 
gously for the notations bri+hep, hep+fee, etc. In Fig. 
[TJ the S and C atoms are distinguished by ball symbols 
of different sizes. The symbol hcp-hcp in Table [TJ indi- 
cates that the S atoms of both, adsorbates 1 and 3, are 
placed exactly at a hep center too (Fig. [TJ;), with analo- 
gous definitions for the symbols bri-bri, top-top etc. The 
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arrangements top+hcp' (Fig. [TJl) and top+hcp differ by 
the dihedral angle of the two S-C bonds, and accordingly 
for other primed structures. The top+fcc+3b configura- 
tion is shown in Fig. [JJ3, where the S-S group forms three 
bonds with the surface. 

From TableQJ at 0.5 ML the adsorption energy for the 
most stable structure is 0.48 eV (top+top, see Fig. [Tg), 
and the S adsorption sites are on the top of Au atoms, 
in keeping with recent experimental results [l(]| . The S-S 
distance is about 3.6 A. Thus at 0.5 ML coverage a dimer 
is formed when methanethiol is adsorbed to the regular 
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Au(lll) surface at low temperature. This dimerization 
was indeed observed recently by STM experiment at 5 
K 0. In Tables HI El and EH 

some configurations at 
1.0 ML coverage are marked by These structures are 
not favored since they involve broken S-C bonds. From 
Table U the adsorption energy of the stable configuration 
at 1.0 ML, with top+hcp as initial structure, is 0.26 eV. 
The optimized structure involves one S atom in the fee 
region but tending toward the atop site; the other S atom 
locates in the fee region, leaning toward the bridge site, in 
agreement with experimental findings 11| . The resulting 
S-S distance is about 3.6 A. This implies methancthiol 
dimer formation also at 1.0 ML coverage. 

The bond between the methancthiol and Au(lll) in- 
volves an ionic contribution due to electron transfer pro- 
ceeding from the surface to the adsorbate. This inter- 
action is well described by the density functional theory 
(DFT) [§j|. Our calculations yield an effective charge 
of -0.2 e, per methanethiol molecule. The binding en- 
ergy of two pure methanethiol molecules, maintaining 
their spatial relation as realized on the surface, results 
as 0.03 cV by DFT [16], 0.02 eV using a coupled cluster 
with single, double and triple excitations (CCSD(T)/6- 
311+G(d,p)) [27j], and the Van der Waals interaction 



between them is 0.06 eV (27j. Being higher than the 



CCSD(T) value by 50%, the DFT prediction for the inter- 
action between two pure molecules is inaccurate, which 
is consistent with recent observations [28|, 129J. However, 
the interaction that leads to dimerization is a surface me- 
diated effect, which strongly exceeds the Van der Waals 
attraction between the dimer constituents [26[. To eval- 
uate the surface mediated interaction between two dimer 
constituents, we study the case of 0.5 ML where the 
preferred configuration is the optimized top+top struc- 
ture (Fig. [Tg). Specifically, we consider the difference 
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3.58) - E 



{top— top) 



(d 



s-s 



5.08), which is interpreted as a reorganization energy, 
namely the energy gain due to dimer formation by re- 
arranging an initially even adsorbate distribution on the 
surface. This energy is 0.24 eV, and thus twelve times 
(0.24 eV/0.02 eV(CCSD(T))) larger than the weak inter- 
action between two pure methanethiol molecules, main- 
taining their spatial relation as adopted on the surface. 
On the other hand, the interaction energy among two 
molecules and Au(lll) is found to be 0.96 eV (top+top). 
The Van der Waals component is thus estimated to con- 
tribute only 6.3% (0.06 eV/0.96 eV) to the total inter- 
action energy. Therefore, the methanethiol dimer forma- 
tion is not caused by the Van der Waals force. Rather, it 
is a surface mediated effect. Since the energy difference 
between the dimer structures at a S - S distance of 3.6 
A0, Q and the reference situation of methanethiol ad- 
sorbates spread evenly over the Au(lll) surface (dss 
around 5.0 A) [lj| is minimally 0.24 eV at all coverage 
levels, for both methancthiol and methylthiolate, omis- 
sion of the relatively small Van der Waals contribution 
does not change our conclusions regarding the dimeriza- 
tion effect. 



B. Methylthiolate adsorbates on the regular 
Au(lll) surface 

Does methylthiolate exhibit the same dimerization be- 
havior as methanethiol? Table |TT] lists the configurations 
of methylthiolate on clean Au(lll) at two levels of cov- 
erage, 0.5 ML and 1.0 ML, with the S-S distances close 
to 2.1 A, 3.6 A and 5.0 A. The maximum adsorption 
energies for 0.5 ML and 1.0 ML, underlined in Table HT1 
are 2.12 eV and 1.99 eV (per molecule). Table HT1 shows 
that at 0.5 ML, the S-S distance is about 3.6 A, which 
is indicative of a dimer (see Fig. Q}i). At 1.0 ML, an 
S-S distance of 5.0 A is found. In the light of this re- 
sult, the absence of dimerization reported in [8j appears 
plausible. The respective measurement involved a cover- 
age level of 1.0 ML, and from our simulation, no dimers 
emerge in this case. The lack of dimerization at 1.0 ML 
for methyltiolate as opposed to methanethiol is rational- 
ized by the much stronger substrate-adsorbate interac- 
tion for the former than for the latter species, as shown 
by the entries for E a d s and ds-Au in Tables HI and HT1 



C. Methylthiolate adsorbates on the defected 
Au(lll) surface 

To examine if substrate defects induce S-S dimer for- 
mation, we allow for a vacancy in the top layer of the 
(3 x 2V3) supercell @ (see Fig. QJ). The location 
of this vacancy is chosen such that it affects as many 
methylthiolate adsorbates as possible. It is sufficient to 
include methylthiolate in this segment of our study since 
methanethiol reduces to methylthiolate when adsorbed 
to the defected Au(lll) surface [HEl. Table [m] com- 
prises the methylthiolate configurations in the presence 
of a vacancy at coverage levels of 0.5 and 1.0 ML. 

With respect to dimer formation, the same observa- 
tions are made as for the regular surface. The dimeriza- 
tion behavior of methylthiolate is therefore not impacted 
by the presence of defects. 



D. Molecular Dynamical Simulation 

In the following, we inspect the effect of tempera- 
ture on methylthiolate dimerization. Two methylthio- 
late molecules are adsorbed on the top layer of one-half 
of the (3 x 2-\/3) supercell, the resulting coverage being 
1.0 ML. Fig. [3 shows the time variation of the S-S dis- 
tance as resulting from MD simulation for a) methylth- 
iolate on the regular Au(lll) at 300K, b) at 400K, c) 
methylthiolate on the imperfect Au(lll) surface at 400K, 
d) methanethiol on the regular Au(lll) surface at 240K. 
At 300K, methylthiolate is still attached to the regular 
Au(lll) surface, the S-S distance oscillates around 5.0 
A, which shows that the adsorbates do not dimerize at 
this temperature. As one raises the temperature to 400K, 
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TABLE II: The geometries and adsorption energies for 
methylthiolate adsorbed on the regular Au(lll) surface at 
0.5ML and 1.0ML coverage. 



TABLE III: The geometries and adsorption energies for 
methylthiolate adsorbed on the Au(lll) surface with a va- 
cancy. 
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however, the S-S distance reduces to 3.5 A. The temper- 
ature increase weakens the molecule-surface bond and 
thus favors intermolecular bonding. This process leads 
to dimer formation. If the temperature is high enough 
and methylthiolate starts to desorb, two molecules join 
to form dimethyl disulfide (CH 3 SSCH 3 ). Thus the ob- 
served dialkyl disulfide emerges during the desorption 
process, and does not exist on the surface @, 0. At 
400K, methylthiolate forms a dimer on regular Au(lll), 
while the S-S distance turns out to be 5.0 A for imperfect 
Au(lll), and no dimerization is found (Fig. [He)). In the 
presence of a vacancy, the interaction between methylth- 
iolate and Au(lll) is stronger than in the regular case 
[1; El, H3| j and interaction between the molecules is ac- 
cordingly weaker. For methanethiol, no dimers exist at 
240 K. 



IV. SUMMARY 

In conclusion, we have demonstrated for the first time 
that at low temperature, the methanethiol S atoms form 
dimers when adsorbed on the regular Au(lll) surface 
pr| . As the temperature rises, these dimers dissociate. 
The S atoms in methylthiolate do not dimerize at 1.0 
ML, while a dimer emerges as most stable at 0.5 ML. 
With increasing temperature the S-S distance of two ad- 
jacent methylthiolate molecules shortens, leading to the 
formation of the experimentally observed dialkyl disul- 
fide. This chemical dimer emerges during the desorption 
process and does not exist on the surface. Prior to this 
study, it was assumed that methanethiol and methylth- 
iolate realize the same pattern upon adsorption on the 
Au(lll) surface. This assumption was based on the 
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FIG. 2: Time variation of the S-S distance in four cases. 



premise that methanethiol deposited on a gold substrate 
turns into methylthiolate. This proposition, however, has 
turned out to be fallacious. As demonstrated by both 
experiment and theory (ill . Il2| , methanethiol persists on 
the clean Au(lll) surface, and this work establishes that 
methanethiol and methylthiolate display distinctly differ- 
ent dimerization patterns. Our results thus resolve the 
long-standing controversy over the correct interpretation 
of the experimental data related to the dimer formation 
by methanethiol and methylthiolate headgroups on the 
Au(lll) surface. 
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